Mechanical instabilities of biological tubes. 
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We study theoretically the shapes of biological tubes affected by various pathologies. When ep- 
ithelial cells grow at an uncontrolled rate, the negative tension produced by their division provokes a 
buckling instability. Several shapes are investigated : varicose, enlarged, sinusoidal or sausage-like, 
all of which are found in pathologies of tracheal, renal tubes or arteries. The final shape depends 
crucially on the mechanical parameters of the tissues : Young modulus, wall-to- lumen ratio, homeo- 
static pressure. We argue that since tissues must be in quasistatic mechanical equilibrium, abnormal 
shapes convey information as to what causes the pathology. We calculate a phase diagram of tubular 
instabilities which could be a helpful guide for investigating the underlying genetic regulation. 
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Tubular structures are found ubiquitously in living or- 
ganisms, from worms to humans and are fundamental 
structures of many organs: they convey liquids, gases or 
cells throughout the body. Typical examples are arteries, 
intestinal tubes or renal excretory canals. The disruption 
of these tubes is at the origin of many pathologies, and al- 
though a considerable amount of research has focused on 
the molecular mechanisms underlying each disease and 
each type of tube [1-5], little attention has been given 
to the biomechanical aspects of tube stability. Here, we 
present a general framework to describe the instabilities 
of cellular tubes; we also discuss the possibility to in- 
fer the underlying causes of pathologies by studying the 
shapes of the abnormally deformed tubes. Indeed, tube 
instabilities are observed in many different organs but al- 
ways seem to fall into a limited number of physiological 
categories : enlarged, sinusoidal, or varicose tubes [2] [6]. 
From a physicist's point of view, these shapes are remi- 
niscent of classical fluid instabilities. Nevertheless, a key 
difference from passive systems is that pathological in- 
stabilities result from the internal activity of the tube. 
Cell division, or active fluid pumping, are often the mo- 
tor of the instability. Thus, a mechanical description of 
the tube shapes could bear a lot of information on the 
underlying mechanisms of tube formation and instability. 

In the past years, progress has been made in the me- 
chanical description of tissues as dividing elastic media [7- 
] . Stresses inside the growing tissues are essential to de- 
termine the final architecture of an organ and conversely, 
there are feedback mechanisms of form on growth. Here, 
we apply these models, which have been usually used 
for bulk tissues, to cylindrical geometries. We consider 
a wide variety of epithelial tubes that all share key fea- 
tures. The tube is composed of a layer of dividing ep- 
ithelial cells, which exert pressure on the surrounding 
medium, and of an elastic basement membrane which 
provides mechanical stability. A softer, visco-elastic con- 
nective tissue, which we consider infinite, is surrounding 
the tube. 



This simple description suggests a competition be- 
tween several forces. The fluid-tube interface has a ten- 
sion and can undergo a Rayleigh-Plateau like instability 
[10], reminiscent of a cylindrical jet breaking into drops. 
However, the tube is elastic, and can resist this defor- 
mation as discussed for passive lipid tubes [11]. Here we 
wish to add active effects: epithelial cells grow (some- 
times at an uncontrolled rate), and the tube can buckle 
under the load induced by growth. This can be modeled 
by a uniform negative surface tension: cell divisions are 
driving an increase of the surface area of the tube. As 
noted in [ ], sufficient mechanical feedbacks can explain 
why growth is spatially uniform. We will make that as- 
sumption here, given the regularity of the patterns we 
wish to explain. Tubes can also dilate because of an ex- 
cess uniform fluid pressure. 




FIG. 1. Left: Sketch of our model for biological tubes. A sin- 
gle layer of epithelial cells rests on an elastic membrane (blue) , 
surrounded by soft connective tissue (not drawn). Right: Dif- 
ferent tube configurations: a) reference, b) varicose/pearling, 
c) enlarged, d) sinusoidal e) sausage-like, b) and e) have the 
same symmetry, but correspond to different instabilities, as 
described in the text. 

We successively investigate the morphologies sketched 
on Fig.l, in that order. The Young modulus E t for bi- 
ological tubes, such as arteries [ ], is of the order of 
10 4 — 10 6 Pa, which is stiff compared to the surrounding 
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tissues. The reference state is a tube of infinite length 
and radius 7*0 . We first discuss a peristaltic perturba- 
tion (varicose, or pearling) which respects the cylindrical 
symmetry. The displacements can therefore be written 
(u r (z),u z (z)) where z is the coordinate along the tube 
and r the radial direction. The forces involved are equiv- 
alent to those derived from the following effective energy 
[13]: 

(1) 

where v is the Poisson ratio, and the strains are defined 
as 

e e = — + -(— ; e z = u z + -{u r +u z ) 2) 
r 2 r 2 

where the prime denotes a derivative with respect to z. 
The first term is the stretching elastic energy of the tube. 
We use here the classical Foppl-Von Karman approxima- 
tion [ ], neglecting u' z 2 , but we keep the nonlinear term 
(^) 2 , if uniform dilations of the tube are non-negligible. 
The second term is the layer bending energy where C is 
the local curvature, and K = 12 fi!l u 2^ the bending mod- 
ulus. A key difference with Ref. [ ] is that we do not 
assume that the outer surface of the cylinder is fixed. The 
case of arteries is complicated by the presence of smooth 
muscles and anisotropy but as we wish to capture the es- 
sential physics, we restrain ourselves to this simple non- 
linear and isotropic elastic theory. Surface and pressure 
forces, which drive the deformation, are derived from an 
effective energy in which the non-equilibrium aspect is 
hidden in the surface tension and excess pressure terms: 



8 a = J (2^ 7 (r + u r (z))y/l + u' r (z) 2 

-7rAP(r + u r (z)) 2 )dz (3) 

where 7 is the effective surface tension, which contains 
a negative contribution due to cell division. As a result, 
7 can be positive or negative. AP is an excess pressure 
inside the tube. We first investigate a deformation of the 
tube to a radius u r {z) = ro[Ro — 1 + Acos(kz)] which 
includes both a uniform dilation to a radius Ro^o and a 
pearling deformation of amplitude Aro <C 7*0. All lengths 
have been normalized by ro (see SI for details). We first 
discuss the effect of negative surface tension, assuming 
that there is a small excess pressure, APro <C \^\(h/ro) 2 . 
Compared to the classical Rayleigh-Plateau instability, 
the wavelength observed in vivo is small, close to the 
value of the tube radius [2, 6, 15]. We propose that this 
pearling instability is due to a buckling of the tube in- 
duced by the homeostatic pressure of the dividing ep- 
ithelial cells. This hypothesis has biological grounds for 
hepatic arteries, since varicoses are associated with Fi- 
bromuscular Dysplasia (FMD) [15-17] or polycystic dis- 
orders in renal canals [18]. 



In the following, we assume that the value of the ten- 
sion is such that \ <C 4?j- <C 1, which seems to be the 

Tq till 

case for reasonable values of the parameters : tubes we 
,2 

will study have ^ ~ 0.01, and as we will discuss below, 

r o 

-gr^ ~ 0.1 . The minimization of the energy with respect 

to Ro leads to Ro = 1 + -j^ — 1. In this limit, there is a 
bifurcation from tubular to pearling states for 

l7l > 7p = Eth-^- (4) 
V3r 

at a wavelength X p = 27rro(jj^a v /^ r )^. The finite 
elasticity of the surrounding medium does not change 
these equations as long as its modulus Es is such that 
Es r o <^ Eth, which is the case here. The Poisson modu- 
lus v only has a weak effect. We have performed a numer- 
ical minimization and Fig. 2 shows a plot of the resulting 
shape of a tube that has undergone a pearling instability. 
For large arteries such as the common hepatic artery, the 
wall to lumen ratio is approximately y- = 0.1, so that 
X p w l.lro- Strikingly, clinical studies [ ] confirm that 
the size of the pattern is proportional to the radius of 
the artery, with the exact same coefficient as above, in 
contradiction to analyses based on the Rayleigh-Plateau 
instability [19]. 

Away from the bifurcation point, the wavelength varies 
only slowly. This variation is compatible with the biolog- 
ical variability, and explains why the same range of wave- 
lengths is always observed. Since the wavelength is of the 
order of the radius, we need not consider instabilities in 
the section plane of the tube, studied in Ref. [ ]. For 
much thinner tubes, such an analysis would be required. 
The negative tension exerted by the cells is linked to the 
homeostatic pressure, 7 « —Phh [ ]. Pressure measure- 
ments for tumors or dysplasia often give values in the 
range of P h = 10 4 Pa [ ]. Taking E t = 10 5 Pa, the crit- 
ical tension j p w 7.10 3 h is perfectly compatible with a 
buckling induced by the dividing cells. 
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FIG. 2. Comparison between the shapes obtained from buck- 
ling theory and in vivo data. A: a common hepatic artery 
displaying varicosities [15]. B: a drosophila tracheal tube dis- 
playing sinusoidal mutation [ ] , A « 5ro . 

Polycystic diseases are more complex. The tubular 
canals of the kidney loose their stability, giving rise to 
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a variety of patterns : string of small pearls, large con- 
volutions, or huge spheric cysts that can invade the en- 
tire organ [6] [24]. Three mains factors [18] have been 
identified : increase in basement membrane compliance, 
uncontrolled division of the epithelial cells forming the 
tube wall, and disorders in ions pumps, causing excess 
water to be pumped in the canals. In our formalism, 
this relates respectively to a decrease in E t , an increase 
in P^, and an increase in AP. The excretory canal of 
C. Elegans is often used as a model organism to study 
these diseases. Although the tube is different in nature 
[25] (a single, elongated cell invaginates to form a tube 
which spans its entire length), it can be described by 
similar equations : there is a bending energy associated 
to membrane deformation, a surrounding elastic medium 
prevents stretching, the tube can grow because of exces- 
sive lipid transport, analogous to negative surface ten- 
sion, and there can exist hydrostatic pressure differences. 
Thus, in order to incorporate the effect of pressure dif- 
ferences, we now consider the opposite limit of enlarged 
tubes, where the pressure dominates over tension effects, 
APro ^> \j\(h/ro) 2 . The equilibrium radius is such that 
the elastic deformation equilibrates the excess pressure. 

The physical picture of pearling is not modified greatly 
by considering the effect of an excess pressure. The 
threshold tension remains of the same order, and the 
characteristic undulation wavelength is slightly reduced 
: X p = X P (AP = 0)(1 + 5APr /3E t h)-^ 4 . 

We then turn to sinusoidal oscillations of the entire 
tube. The classical Euler buckling instability of an elastic 
column under compression occurs at a wavelength which 
is the system size, in the absence of a surrounding elastic 
medium. However, in an elastic medium, buckling occurs 
at a finite wavelength. We have shown recently that this 
instability can be invoked to describe the morphology of 
the intestine [26]. 

We call E t and E s respectively the Young moduli of the 
tube and surrounding elastic medium, keeping in mind 
that E t » E s . We denote the deformation of the tube 
by w(z,t) = Bcos(kz). The curvature modulus of a hol- 
low elastic tube is K = TrE^k, the force exerted by 
the growing tube is / = 27T7ro, and the elastic force ex- 
erted by the surrounding medium is proportional to the 
deformation with a coefficient a. The threshold force of 
the sinusoidal instability is / = 2\JKol^ thus the critical 
tension for buckling is 

| 7 | > 7s « 0.8ViJiVE s E t (6) 

The wavelength of the instability is X s = 27rro ( Jr 1 ) ^ . 
A rough estimate leads to A c « 5 — 10 ro, in agreement 
with observations of degenerate excretory canals in C. 
Elegans [ ], or of the tracheal system in drosophila ([ ] 
and Fig. 2). 



We now compare sinusoidal deformations, which oc- 
cur for cell tensions larger than j s given by Eq.(6), 
with pearling deformations with a threshold j p given by 
Eq. (4) . The preferred morphology essentially depends on 

E 

the value of the parameter r — ^p-. If the substrate is 
very soft or the tube very hard, the sinuous morphol- 
ogy is expected, since in that case, it requires a smaller 
tension. Conversely, if ^ is very small, the pearling in- 
stability is expected. Physically, this is due to the fact 
that a larger tube is harder to bend than a thicker tube 
: its curvature modulus varies as the cube of the radius 
and as the first power of thickness. For a typical tube, 
^ =0.1, and E t = 10 5 Pa, the critical substrate rigidity 
separating both morphologies is E s = 10 2 Pa. This anal- 
ysis is important since a loss of tubular rigidity is one 
of the causes of cystic diseases. The situation is much 
more complex if the substrate is modeled as viscoelastic. 
Over long timescales, sinuous deformations are always 
expected, but can be kinetically constrained. A detailed 
discussion goes beyond the scope of this paper [27]. 

In the region where both instabilities can be observed, 
it is necessary to compare the energies of the two con- 
figurations. We plot in thick green the transition line 
between the two instabilities, in the vicinity of the criti- 
cal point. 




FIG. 3. Stability diagram of all possible configurations of a 
tube, depending on the relevant parameters. The vertical axis 
is the ratio of the moduli of the surrounding tissue and of the 
tube. The horizontal axis is the ratio of the tension, negative 
or positive, and the elasticity of the tube. Lines in thin red 
represent the transition to an instability. The thick green 
line separates stability domains of the varicose and sinuous 
patterns. 

The results are displayed on a stability diagram in 
Fig. 3. The drosophila tracheal system would be a good 
candidate for putting this theory to the test. It has been 
widely studied as a model of tube size regulation and 
the deletion of several genes has been identified to cause 
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pathologies. Specifically, mutants lacking the Sinuous, 
or Varicose genes respectively display sinuous [ ], and 
pearling state morphologies [ ]. Although we do not 
have a good understanding of their roles, it is known that 
they are required in cell-cell junctions (associated to the 
rigidity of the tube) and that they have a role in control- 
ling growth. We now show how our phase diagram could 
guide further experiments. The modulus of the normal 
tracheal wall was inferred from experiments to be around 
20kPa[ ], and ^ w 0.2. Although tracheal connective 
tissue has not been probed, usual orders of magnitude are 
lOOPa. In normal organisms, the relevant ratio therefore 
is = 0.6 < 1, so sinusoidal instabilities are expected. 
Nevertheless, disturbing junctions lowers the rigidity of 
the tube, and a three-fold decrease in E t would cause the 
pearling instability to be favored. We therefore suggest 
that more emphasis should be put on measuring each 
gene influence on the mechanical properties of cellular 
tubes. 

A last situation to address is a tube in the absence 
of any dysplasia of the cell wall. The surface tension is 
then positive. In this regime, we expect either a sta- 
ble tube, or a classical Rayleigh-Plateau instability of an 
elastic cylinder. This is slightly different from the situa- 
tion encountered in membrane tubes [ 1], since the fluid 
cylinder is surrounded by an elastic shell. 

We study a perturbation u r (z) = ro[(Ro — 1) + 
Acos(kz)]. The situation that we wish to discuss does 
not involve any cell growth, and the instability thus 
occurs on a time scale small enough volume conserva- 
tion is a reasonable assumption. This imposes that 
Rq — 1 = — 4>~- Minimizing the energy with respect to 
A (SI) yields a finite value of A only for a tension larger 

i — 2 

than the critical value j r = E t h j_^ 2 . Below this critical 
value of the tension, the cylinder is stable. The instabil- 
ity first occurs for the longest wavelength mode q — » 0. 
Determining the exact most unstable wavelength requires 
a non-linear analysis, which is beyond the scope of this 
paper. 

Nevertheless, if the size of the perturbation approaches 
the radius of the tube, disconnected pearls are obtained 
and the most unstable wavelength can be found by study- 
ing the dynamics of the perturbation u z . The continu- 
ity equation and a Poiseuille approximation for the flow 
lead at lowest order to: = r^/ (lQr])d 2 P where r] 

the fluid viscosity. The hydrostatic pressure P is calcu- 
lated from the force-balance equation on the membrane 
P oc (jA(-q 2 + 1) - E t h(l - v 2 )A) cos(qz). The most 

unstable wavenumber is then q 2 . = ^- Eth ^~ v ) m This 
is an explanation similar to [19] for the sausage-like pat- 
tern observed on some arteries. Indeed it has been shown 
[31] that a very soft elastic cylinder can undergo a classi- 
cal Rayleigh-Plateau instability. Testing the case of hol- 
low elastic tubes in the same fashion could be an ex- 
perimental confirmation to our theory. As expected, if 



j/E t h ^> 1, the most unstable wavevector converges to 
the wavevector of the classical Rayleigh-Plateau instabil- 
ity q R p = ^ [10]. 

The main result of this paper is the stability diagram, 
sketched on Fig. 3 presenting the stabilities of most mor- 
phologies found in cellular tubes, depending of their me- 
chanical properties. Buckling theory is in good quanti- 
tative agreement with the in vivo data. We have here 
deliberately used simple models for the mechanics of the 
tubes, to emphasize the generality of our results, and the 
relevance of a few dimensionless ratios to predict the ob- 
served instabilities. Most importantly, we want to stress 
the fact that tissues' shapes could turn out to be a mean- 
ingful criterion to detect the causes of specific patholo- 
gies. For example, in the case of excretory kidney canals, 
dysplasia leads to a different instability than fluid accu- 
mulation. More work is clearly needed in each specific 
case, but we do not expect the basic physics presented in 
this paper to be modified. 
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